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Abstract 

We study the behaviour of the flux tube thickness in the vicinity of 
the deconfinement transition. We show, using effective string methods, 
that in this regime the square width increases hnearly and not loga- 
rithmically with the interquark distance. The amplitude of this linear 
growth is an increasing function of the temperature and diverges as the 
deconfinement transition is approached from below. These predictions 
are in good agreement with a set of simulations performed in the 3d 
gauge Ising model. 



1 Introduction 



One of the most intriguing features of the confining regime of Lattice Gauge 
Theories (LGTs) is the logarithmic increase of the square width w'^{R) of 
the flux tube as a function of the interquark distance R [Jj. This effect, 
which is commonly referred to as the "delocalization" of the flux tube is 
together with the well known Liischer term one of the most important and 
stringent predictions of the effective string approach to LGTs. 

This effect was discussed for the first time many years ago by Liischer, 
Miinster and Weisz in [T] but it required several years of efforts before it 
could be observed in lattice simulations. The first numerical results were 
obtained in abelian models [2H5] and only recently they were extended also 
to non abelian LGTs [6l[7] (see also pj for some early attempt). In all these 
tests a good agreement with the theoretical predictions was found. 

A question which naturally arises in this context is which is the fate 
of the flux tube (and in particular the behaviour of its width) as the de- 
conflnement transition is approached from below. It is important to stress 
that delocalization and deconfinement are two deeply different conditions 
of the flux tube. The deconflnement transition is characterized by the van- 
ishing of the string tension cr(T). In pure gauge theories it is associated 
to the breaking of the center symmetry of the gauge group and its natural 
order parameter is the Polyakov loop. At the deconfinement point the flux 
tube vanishes. The delocalization of the flux tube instead coincides with the 
roughening transition. In all the physically interesting models this transition 
occurs for value of the gauge coupling (3 well below the fixed point where 
the continuum limit can be taken. Thus in the continuum limit (for any 
temperature below the deconfinement transition) the theory is always in the 
rough phase and thus the fiux tube is always delocalized. Delocalization is a 
typical quantum effect. It is a consequence of the Mermin- Wagner theorem 
which imposes the restoring (in the continuum limit) of the translational 
symmetry for the fiuctuations of the flux tube in the transverse directions. 
Intutively it amounts to say that we cannot flx deterministically the trajec- 
tory of the flux tube but may only describe it as a probabilty distribution. 
It is important to stress that, even if delocalized, the flux tube fully keeps 
its conflning function. The quantum fluctuations which drive the delocaliza- 
tion also influence the conflning potential (as the presence of Liischer term 
indicates) but do not destroy it. 

While the behaviour of the string tension cr(T) as the deconfinement 
temperature Tc is approached from below is rather well understood much 
less is known on the behaviour of the flux tube thickness in this regime. 
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This is a rather important issue from a physical point of view since the 
interplay between delocalization and deconfinement could strongly influence 
the transition from hadrons to free quarks as Tc is approached. 

This problem can be addressed within the framework of the effective 
string approach by performing a modular transformation of the low temper- 
ature result. This was done in [9| in the case of the free bosonic approxima- 
tion (i.e. the first order in the perturbative expansion of the Nambu-Goto 
effective string). In this free bosonic limit one can show that when peri- 
odic boundary conditions are imposed in the time direction (i.e. when finite 
temperature regularization is imposed) the large R behaviour of the square 
width changes completely and becomes linear instead of logarithmic. This 
behaviour holds in principle for any temperature T, but as T decreases it 
requires larger and larger values of R to be observed. Similarly it is possible 
to show that for any fixed value of R the square width smoothly converges 
toward the expected logarithmic behaviour as T decreases. The threshold 
between the two behaviours is i? ~ 1/T. The coefficient of the linear de- 
pendence at high T is itself temperature dependent. It increases linearly 
with the temperature and remains finite as the deconfinement transition is 
approached. 

In [9j these predictions were compared with a set of high precision simu- 
lations in the 3d gauge Ising model, with results which turned out to be only 
in partial agreement with the effective string picture. For all the temper- 
atures studied in [9] vP'{R) was indeed, (with very good confidence levels) 
a linearly increasing function of R. However the coefficient of this linear 
behaviour was in general larger than that predicted by the effective string 
(except for the smallest temperature values) and, what is more important, 
it seemed to diverge as the deconfinement point was approached. 

It is clear that in order to understand this discrepancy one should go be- 
yond the first order approximation in the effective string calculation. As a 
matter of fact in this finite temperature description, the perturbative expan- 
sion in powers of l/a^Lp', is actually an expansion in T'^ jT'^ (see below for 
further details on this point) thus it is somehow obvious that higher orders 
terms become important if one wants to study the behaviour in the vicinity 
of the deconfinement point. What is less obvious is how relevant are for the 
fiux tube width these higher order contributions. Indeed it seems that the 
particular problem that we are addressing here, the high temperature be- 
haviour of the fiux tube thickness, is one of the best numerical laboratories 
to test higher order contributions in the effective string description. In fact 
the comparison with the 3d Ising model [9] shows that these higher order 
corrections may become very large (much larger than the first order term) 
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in the vicinity of the deconfinement transition. 

Higher order effective string corrections have attracted much interest in 
these last years both as a tool to better fit the numerical results and as 
a way to improve our understanding of the interplay between the effective 
string approach and standard string theory. However no result had been 
obtained until very recently on higher order corrections to the flux tube 
width. Most of the papers concentrated instead on higher order corrections 
to the interquark potential. Several important results were obtained in this 
context. It was shown that not only the leading (Liischer) term but also the 
first subleading (quartic) correction is universal |10H12] . thus remarkably 
enhancing the predicting power of the effective model. In a recent paper 
universality has been extended also to the subsequent term (sixth order 
term in the potential) |13j at least for the (2+1) dimensional case and for 
the torus and cylinder topologies. Finally, in the case of the Nambu-Goto 
string, the all order expansion was obtained in |10pi4j. 

The situation changed recently thanks to [6] where an explicit expres- 
sion for the next to leading (quartic) correction to the flux tube width was 
obtained in the case of the cylinder topology. This is exactly the topology 
which is needed to describe the behaviour of the flux tube at finite temper- 
ature. 

The aim of the present paper is to use the result of [6j to address the 
discrepancy between theoretical predictions and numerical results described 
above. As we shall see, keeping into account this correction the right T 
dependence of the flux tube thickness is recovered. This, together with a 
parallel result obtained using a dimensionally reduced effective model for 
the high temperature behaviour of lattice gauge theories, will allow us to 
guess further terms in the expansion and to propose the following conjecture 
for the behaviour for the large R limit of the flux tube square width up to 
the deconflnement transition: 



where cr(T) denotes the temperature dependent string tension. 

If, in addition, one also assumes a Nambu-Goto effective action for the 
flux tube (assumption which we know to be wrong, but in several cases gives 
a very good approximation, at least for the first few orders in T) then the 
T dependence of cr(T) can be predicted exactly and turns out to be 




(1) 
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Thus allowing a complete description of the T dependence of the flux tube 
thickness in the limit of large interquark separation. 



2 Effective string prediction for the flux tube thick- 
ness at flnite temperature 

2.1 Definition of the flux tube thickness 

In a finite temperature setting the lattice operator which is used to evaluate 
the flux through a plaquette p of the lattice is: 



where P, P' are two Polyakov loops separated by R lattice spacings and Up is 
the operator associated with the plaquette p. Different possible orientations 
of the plaquette p measure different components of the flux. In the following 
we shall neglect this dependence which plays no role in our analysis. The 
only information that we need is the position of the plaquette. Let us define 



where h denotes the displacement of p from the P P' plane. In each trans- 
verse direction, the flux density shows a gaussian like shape (see for instance 
Fig. 2 of [2] ) . The width of this gaussian w is the quantity which is usually 
denoted as "flux tube thickness" : 



This quantity depends on the number of transverse dimensions and on 
the bare gauge coupling /3 (or, equivalently, on the lattice spacing a). Once 
/3 is fixed the only remaining dependences are on the interquark distance R 
and on the lattice size in the compactified timelike direction L, i.e. on the 
inverse temperature of the model. By tuning L we can thus study the flux 
tube thickness near the deconfinement transition. 




(3) 



{cPip;P,P')) = (^cP{h-R,L)) 




(4) 
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2.2 Effective string prediction: general setting 

In the effective string framework the square of flux tube width is given by: 

/ [Vh] h\t,x)h\t,x)e-^^^^ 

w'^{x; R, L) = (5) 

J [Vh]e-^^^^ 

where the sum is intended over ah the surfaces bordered by the two Polyakov 
loops and S[h] is the effective action. As mentioned above, up to the second 
order in which we are interested here, this action is universal and coincides 
with the expansion to the fourth order in h of the Nambu-Goto action. In 
the following we shall concentrate in the (2+1) dimensional case (i.e. we 
shall deal with only one transverse direction) and eliminate the index of the 
field h. In this case the Nambu Goto action has a particularly simple form 

rL/2 i-R/2 

S[h] =ao dT dq^Jl + {drhY + (5,/i)2 . (6) 

J-L/2 J~R/2 

The correlator eq.dSJ is singular and must be regularized. The standard 
choice is the point splitting regularization. It is easy to see looking at eq.([6]) 
that, in the large R/L limit in which we are interested there is a natural ex- 
pansion parameter in the evaluation of eq.([5]) which is l/{aoLp'). The leading 
order in this expansion is simply given by the free bosonic approximation 
to the Nambu-Goto action while the next to leading correction will include 
the fourth order terms of the action. Let us address these two contributions 
in more detail. 



2.3 Effective string prediction: leading order 

The leading order effective string prediction for square of the flux tube thick- 
ness wf^ in the finite temperature (i.e. cylinder) geometry was evaluated for 
the first time (in a rather implicit form) in [2]. It was then reobtained in a 
slightly more explicit form in ^ using the method of images. We shall use 
as our starting point this last expression: 

^o^Fo = + i 1^2(0)/^; (0)1 (7) 

where ctq denotes the zero temperature string tension and L is the length of 
the cylinder in the compactified direction (i.e. the inverse temperature) and 
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the square width is evaluated at the midpoint between the two Polyakov 
loops. 

It is easy to see using standard relations among 9 and rj functions that 
this expression coincides with the one reported in 

In the large R and high temperature limit, i.e. for R» L this expres- 
sion becomes: 

which shows, as anticipated, that the square width should increase lin- 
early with R with a coefficient which increases linearly with the tem- 
perature and stays finite at the deconfinement point. In fact the zero tem- 
perature string tension do (which can be extracted from the Wilson loop or 
from Polyakov loop correlators at very low temperature) does not depend 
on the finite temperature T but only on the coupling constant (3. 

2.4 Effective string prediction: next to leading correction 

The next to leading correction was recently evaluated in |6].We report here 
for completeness the result: 



E2{t) - ^E2{2t) 

= 48 ' 

where E2{t) is the first Eisenstein function: 

oo „ 
n=l 



r = 1^ denotes the modulus of the cylinder and 



q — Q^TTIT _ ^--kL/r 

The dominant term in the R » L limit (i.e. r — )• 0) turns out to be 
again linear in R: 
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(10) 



24(72^3 



+ ••• 



This result is rather non trivial. Looking at the small r expansion of the 
Eisenstein function one would in principle expect a term proportional to B? 
which however disappears in eq.® as a consequence of a set of non trivial 
cancellations among Eisenstein functions. This is an important consistency 
check in view of the fact that the simulations reported in [9j found no evi- 
dence of terms proportional to B? and found instead a remarkably precise 
linear increase of up' with R. 

Combining together the leading and subleading corrections we end up 
with the following expression: 



which, assuming the Nambu-Goto value for the critical temperature = 
Scto/tt (see below for a derivation of this result), can be rewritten as an 
expansion in T/T^ 



We compare in tab{T]and fig{T] (dashed line) this result with the values of 
the numerical simulations reported in [9j. We see a clear improvement with 
respect to the leading contribution alone (continuous line), which however 
is not enough to fill the gap with the numerical data. At the same time 
it is easy to see that, as anticipated, for values of L in the vicinity of the 
deconfinement point the subleading term becomes of the same order of the 
leading one. This shows that in this limit the expansion in powers of l/uo-L^ 
converges too slowly and that it would be very important to have some 
educated guess on how to resum the expansion. We shall show in the next 
section that this guess can be obtained with a completely different approach, 
i.e. looking at the effective theory obtained integrating out the spacelike 
degrees of freedom following the well known Svetitsky-Yaffe proposal (15) . 

2.5 Dimensional reduction and the Svetitsky Yaffe approach. 

In the vicinity of the deconfinement transition the physics of a (d+l) LGT 
can be described using an effective model in which the spacelike links are 
integrated out and the only remaining degrees of freedom are the Polyakov 
loops. This essentially amounts to integrate out all the degrees of freedom 






(12) 
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of the original LGT up to the scale L = 1/T and in fact the theory ob- 
tained in this way effectively behaves as a model in one dimension less than 
the original LGT. It is easy to see that this procedure gives as effective de- 
scription a d dimensional spin model (the Polyakov loops playing the role 
of spins) with a global symmetry represented by the center C of the gauge 
group of the original LGT. This approach was pioneered by Svetitsky and 
Yaffe p!5| who were able to predict in this way the critical behaviour of the 
LGT at the deconfinemenet point under the assumption that both the de- 
confinement transition of the original LGT and the magnetization transition 
of the effective spin model are continuous. However its utility goes beyond 
the prediction of the critical indices. In the last few years it has been shown 
that this effective description can be used to predict the behaviour of various 
quantities in the vicinity of the deconfinement transition as far as the scales 
involved in these observables (for instance the interquark distance R) are 
larger than L = 1/T which sets the scale of the effective theory (for a review 
of most recent results see for instance [E]). This is exactly the R » L 
limit in which we are interested in this paper. 

The simplest examples of this effective mapping are the (2+1) SU{2) 
LGT and the (2+1) Ising gauge model which have the same center Z2 and 
are thus both mapped into the 2d spin Ising model. In this case, since 
both LGTs have a second order deconfinement phase transition we can also 
predict, following [15], that these transitions are in the same universality 
class of the 2d Ising model. 

Let us address in more detail these examples. It is easy to see that the 
confining phase of the LGT is mapped into the high temperature phase of 
the spin model and that, as we mentioned before, the Polyakov loops of the 
LGT are mapped into the spins of the 2d Ising model. This correspondence 
should be understood in the renormalization group sense i.e. the Polyakov 
loop operator is actually mapped into a linear combination of all the (C- 
odd) operators of the spin model. In the 2d Ising case this means that 
we have a combination involving the whole conformal family of the spin 
operator. For R large enough this combination is dominated by the relevant 
operators which in the Ising case is only the spin operator. In a similar way 
it is possible to show that the plaquette operator is mapped into the most 
general combination of the energy and the identity conformal families. This 
allows us to construct the analogous of the operator which measures the 
flux tube thickness which turns out to be a suitable combination of three 
point correlators of the spin and energy operators (see [T7] for a detailed 
discussion of this mapping). In the particular case of the 2d Ising model 
these correlators can be evaluated exactly leading to the following expression 
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for the "flux" distribution [IT] 



27rR e-™VV+«' ^ ^ 

^'"■''> - 4y^ + B.^ A-„(,„Ji) ■ (1^' 

where y denotes the transverse direction, Kq is the modifled Bessel function 
of order 0, m is the mass of the 2d Ising model and a large mR limit is 
assumed. 

From this flux distribution it is easy to extract the square of the flux 
tube width as the ratio 

w\R) = ^-^ / J (14) 

^ ' r^dyP{R,y) ^ ' 

which, setting x = 2y/R amounts to evaluate 



D2 x'^ -2mr^l+x^ 

v^\R) (15) 

4 TOO 7 p — 2mrV l + a:'^ 

These integrals can be evaluated asymptotically in the large mR limit (see |17pl ) 
leading to the following result: 

w'^(R) ~ \ — + .... (16) 
4 m 

where the dots stay for terms constant or proportional to negative powers 
of R. 

The last step in order to compare this result with eq. (|lip is to give a 
meaning to the Ising mass m in terms of LGT quantities. 

This can be easily accomplished if we recall that the mass is the inverse 
of the correlation length of the model and can be obtained from the large R 
limit of the 2 point function of the model as follows: 

Urn {a{0,0)a{0,R)) Ko{mR). (17) 

According to the mapping discussed above this correlator is the 2d limit 
of the expectation value of two Polyakov loops at distance R. In order to be 
consistent with the assumptions we made to evaluate the flux tube width in 
the previous section, we must evaluate this correlator in the framework of 



^Notice, to avoid confusion, that in 17 we used the variable r = R/2 and that we 
evaluated the unnormalized square width, i.e. only the numerator of ea. (|14|) 
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the Nambu-Goto effective action. Tliis result was obtained a few years ago 
by Liisclier and Weisz [10] using a duality transformation and tlien derived 
in the covariant formalism in [l3| . In d = 2 + 1 dimensions one finds a tower 
of Kq Bessel functions: 

oo 

(P(0, 0)P(0, R)) = CnKoiEnR). (18) 

ra=0 

(see eq. (3.2) of [TO]). 

where c„ are constants (whose exact expression is irrelevant for our anal- 
ysis, the only important point is that they do not contain any R dependence 
in (2+1) dimensions, contrary to the higher dimensional cases) and En are 
the closed string energy levels: 



i^. = .oL{l + ^[-^(d-2)+n]} 



1/2 

(19) 



(see eq. (C5) of [TO]). 

It is easy to see that in the large R limit only the lowest state (n = 0) 
survives and, in full agreement with our dimensional reduction mapping, we 
end up with a single Kq function: 

lim (P(0,0)P(0,i?)) ~Ko(^oi?). (20) 

R—>oo 

comparing with eq. p7p we immediately rccognizG 771 — Eq i.e. 
This result can be rewritten as 

m = a{T)L (22) 

where ct(T) is given by eq.© which wew report here: 

2^2 \ 1/2 



a{T) = (^1 - . (23) 

where the critical temperature T"^ = 3<to/71" is the Nambu-Goto prediction for 
the deconfinement temperature and is given by the solution of the equation 
(t(T) = 0. Eq. (j23p is the prediction of the Nambu-Goto effective action for 
the temperature dependence of the string tension and was discussed several 
years ago by Olesen in [ISJ. 
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Plugging this result in eq. ()16p we find for the coefficient of the linear 
term: 



which may be expanded in powers of as follows 

+ + + (25) 

or equivalently 



2 „ 3 T ST^ 



The first two orders exactly coincide with those of eq. ()lip . This repre- 
sents a remarkable consistency check of both the Svetitsky-Yaffe approach 
and the second order perturbative calculation of [6] and leads to conjecture 
that eq. (j24p could be the resummation to all orders in l/{aoL'^) of the lin- 
ear term of the square width in the framework of the Nambu-Goto effective 
action. 

It is important to stress, as a final comment on this result, that even if 
it was obtained assuming for ct(T) the functional form predicted by the 
Nambu-Goto action: eq.([2]) its validity goes well beyond this particular 
model. In fact the calculations which lead to eq. (|lip and to the first two 
terms of eq. (|26p only involve contributions up to the next to leading order in 
the effective string perturbative expansion which were proved to be universal 
in [lOHni. 



2.6 Comparison with the Ising model. 

In figdland tab [T] we compare eq. p4p with the results of a set of simulations 
performed in the 3d gauge Ising model and reported in [9j. In these simula- 
tions the coupling of the Ising model was fixed at /3 = 0.75180 for which the 
deconfinement transition is known to occur at exactly 8 lattice spacings [19] . 
Moreover for this value of /3 the zero temperature string tension is known 
with very high precision to be (Tq = 0.0105255(11) pO] thus excluding a 
possible source of systematic errors. The model was studied for the values 
of L reported in the first column of tab(T] which correspond to temperatures 
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Figure 1: Plot of k{L) as a Junction of = X' '^^^ continuous line is the pre- 
diction for k{L) according to eq.^. The dashed line according to eg. ill]) . The 
dashed dotted line to eg. {2^^ and the dotted line to h25]) . In all the four cases the 
curves were obtained only using the terms explcitly written in the corresponding 
equations and neglecting higher order corrections ( denoted with the dots in the var- 
ious equations). The points are the results of the simulations in the 3d gauge Ising 
model 
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L 


k{L) 


eq.® 


eq-dH]) 


eq.dMD 


eq.jSS]) 
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6.19(10) 


2.639 


4.260 


NAN 


5.753 


10 


4.90(4) 


2.375 


3.557 


33.307 


4.438 


11 


3.85(4) 


2.159 


3.047 


5.121 


3.594 


12 


3.14(4) 


1.979 


2.663 


3.560 


3.017 


14 


2.33(3) 


1.697 


2.127 


2.418 


2.291 


16 


1.84(3) 


1.484 


1.773 


1.899 


1.857 



Table 1: Values of the coefficient k(L) of eg. |^7| ) for various values of L. In the 
second column we list the results of the simulations. In the following columns we 
report the prediction for k{L) assuuming eg. 07]) {24^ ^25 ]) and neglecting higher 
order corrections denoted by the dots in the eguations. 

ranging from T = Tc/2 to T = 8Tc/9. For each temperature w'^ was then 
evaluated for several values of R and fitted with 

w^{R,L) = k{L)R + c{L) (27) 

For all temperatures the reduced of the fits was very good. The values of 
k{L) obtained in this way are reported in the second column of tab{TJ Look- 
ing at figHland tablUwe see that the for all the temperatures eq. (fM|l predicts 
values of k{L) larger than those found in the simulations and diverges for 
T ~ 4Tc/5. This is not surprising. Indeed it is by now well known that the 
Nambu-Goto effective string is not a good description of the 3d gauge Ising 
model. In fact it predicts a deconfinement temperature which is ~ 0.8 of 
the observed one (in agreement with what we find here) and predicts too 
large corrections in the interquark potential. This disagreement was dis- 
cussed in the past years in a set of high precision numerical tests both in the 
high T [21\ (but still in the confining phase) and in the low T [22] regimes 
of the interquark potential and more recently looking at the interface free 
energy of the 3d Ising spin model [20]. This last observable is related by 
duality to the interquark potential and represents an important cross test 
because it involves different (i.e. toroidal) boundary conditions. In all these 
tests the interquark potential was rather well described by the truncation 
at the second order of the Nambu-Goto expansion in powers of l/a^L'^ and 
definitly incompatible with any further higher order contribution. Remark- 
ably enough we see the same pattern also in this case. We report in fig.l (the 
dotted line) and tab.l (last column) the values of k(L) obtained truncating 
eq. ()24p at the second order (i.e. eq. (j25p l. As shown in fig.l the measured 
values lie slightly above this curve. This is indeed a puzzling result in view 
of the recent claim [13] on the universality of the sixth order term in the 
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effective string expansion for the interquark potential. This sixth order term 
seems instead to be very different from the Nambu-Goto one (and compat- 
ible this zero) in the case of the 3d Ising model. We have for the moment 
no good answer to this puzzle. 

3 Conclusions 

Let us summarize our main results and add a few concluding comments: 

1) As anticipated in the introduction our main result is that approach- 

ing the deconfinement transition has a twofold effect on the square 
width of the flux tube. First, the logarithmic dependence on the in- 
terquark distance is changed into a linear one. Second, (see eq.([T2|)) 
the amplitude of this linear growth increases more than linearly with 
the temperature as the transition point is approached. 

2) We conjecture that the behaviour described by eq. (jl2p is only the first 

term of an expansion in powers of T^/T^ and that the square width w'^ 
should be proportional to the inverse of the finite temperature string 
tension (see eq. (|24p ) If we isolate in eq. ()24p the term linear in R we 
find an universal relation which should hold for any LGT in the limit 
of large interquark distances 

a{T)w'^{R) = ^RT (28) 

Moreover, if we assume that the effective action is of the Nambu-Goto 
type, then we may also predict the explicit expression for cr{T) which 
is given by eq.©- 

3) However it is well know that for many models (and in particular for 

the 3d gauge Ising model) the Nambu-Goto picture is not correct and 
that eq.(l2|) is a rather poor description of the actual finite temperature 
behaviour of the string tension cr(T). In the particular case of the 3d 
gauge Ising model we find a good agreement between predictions and 
numerical data if we truncate at the second order the perturbative 
expansion in powers of l/a^L'^ in agreement with what already found 
when looking at the interquark potential. 

4) On the contrary in the case of (2-1-1) SU{N) gauge models (and in 

particular for A'^ = 3) the Nambu-Goto string gives a much better 
description of the finite temperature behaviour of the theory (see |23j ) 
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and accordingly we expect that eq. ()24p should well describe the finite 
temperature behaviour of uP . It would be very interesting to extend 
to these cases the present analysis, using for instance the methods of 
[6] or of [7] and to test on a quantitative basis our prediction eq.([T|). 

5) For LGTs with a first order deconfinement transition (as it happens for 
instance for > 4 in d = 2 + 1 [23] and for SU{2,) in d = 3 + 1), there 
is no reason to expect a divergent behaviour of the coefficient k{L) . 
It would be very interesting to see if eq.([T]) still holds in these cases. 
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